Gaussian flexibility with Fourier accuracy: the periodic von Neumann basis set 
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We propose a new method for solving quantum mechanical problems, which combines the flexibil- 
ity of Gaussian basis set methods with the numerical accuracy of the Fourier method. The method 
is based on the incorporation of periodic boundary conditions into the von Neumann basis of phase 
space Gaussians [F. Dimler et al., New J. Phys. 11, 105052 (2009)]. In this paper we focus on the 
Time-independent Schrodinger Equation and show results for the harmonic, Morse and Coulomb 
potentials that demonstrate that the periodic von Neumann method or pvN is significantly more 
accurate than the usual vN method. Formally, we are able to show an exact equivalence between 
the pvN and the Fourier Grid Hamiltonian (FGH) methods. Moreover, due to the locality of the 
pvN functions we are able to remove Gaussian basis functions without loss of accuracy, and obtain 
significantly better efficiency than that of the FGH. We show that in the classical limit the method 
has the remarkable efficiency of 1 basis function per 1 eigenstate. 

PACS numbers: 2.70.Hm, 2.70.Jn, 3.65.Fd 82.20.Wt 



The formal framework for quantum mechanics is an 
infinite dimensional Hilbert space. In any numerical cal- 
culation, however, a wave function is represented in a 
finite dimensional basis set and therefore the choice of 
basis set determines the accuracy. The optimal basis set 
should combine accuracy and flexibility, allowing a small 
number of basis functions to represent the wave functions 
even in the presence of complex boundary conditions and 
geometry. Unfortunately, these two criteria — accuracy 
and efficiency — are usually in conflict, and globally ac- 
curate methods [iHl] lack the flexibility of local methods 
d, [B| . For example, in the pseudospectral Fourier grid 
method the wave function is represented by its values on 
a finite number of evenly spaced grid points. Due to the 
Nyquist sampling theorem, this allows for an exact repre- 
sentation of the wave provided the wavefunction is band 
limited with finite support0-@|. However, the non-local 
form of the basis functions in momentum space leads to 
limited efficiency. On the other hand, in the von Neu- 
mann basis set [3, [lo| each basis function is localized on 
a unit cell of size h in phase space. However, despite the 
formal completeness of the vN basis set[ll|, attempts to 
utilize this basis in quantum numerical calculations have 
been plagued with numerical errors [iT 
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The purpose of this paper is to establish a precise 
mathematical formalism for the von Neumann basis on 
a truncated phase space. This allows us to demonstrate 
an exact equivalence with the Fourier method while re- 
taining the flexibility of a Gaussian basis. It puts on a 
rigorous basis the seminal work of Dimler et al. [l4| who 
used a vN basis with periodic boundary conditions, al- 
though in our formalism the periodicity of the vN basis 
appears only implicitly. 

The von Neumann basis set [l(| is a subset of the "co- 
herent states" of the form: 
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where n and I are integers. Each basis function is a Gaus- 
sian centered at (na, — ) in phase space. The parameter 
a = 4 s - controls the FWHM of each Gaussian in x and p 
space. Taking Ax — a, Ap = h/a as the spacing between 
neighboring Gaussians in x and p space respectively, we 
note that AxAp — h so we have exactly one basis func- 
tion per unit cell in phase space. As shown in[ll| this 
implies completeness in the Hilbert space. 

The "complete" vN basis, where n and / run over all 
integers, spans the infinite Hilbert space. In any numeri- 
cal calculation, however, n and I take on a finite number 
of values, producing N Gaussian basis functions {gi(x)}, 
i = 1...N. Since the size of one vN unit cell is h, the area 
of the truncated vN lattice is given by 5 vN = Nh. 

In the pseudospectral Fourier method, a function ip(x) 
that is periodic in L and band limited in K = f can be 

written in the following form: tp(x) = Y^ n =i VK' T «)^™( x )i 
where x n — S x (n — 1), and S x = ^ = jr. The basis 
functions {6 n (x)} are given by [TBI ]: 
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which can be shown to be sine functions that are periodic 
on the domain [0,L]. The set {9i(x)} i = 1, ..,N spans 
a rectangular shape in phase space with area of S FGR = 
2LP = 2L^ = Nh. Thus N unit cells in the vN lattice 
and N grid points in the Fourier method cover the same 
rectangle with an area in phase space of: 



5 vN = gFGH = Nh 



(3) 



(Fig. [T]). This suggests that N vN basis functions con- 
fined to this area will be equivalent to the Fourier basis 
set. Unfortunately, the attempt to use N Gaussians as 
a basis set for the area in eq.([3]) (Fig. [1]) is unsuccessful, 
a consequence of the Gaussians on the edges protruding 
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from the truncated space. However, by combining the 
Gaussian and the Fourier basis functions we can generate 
a "Gaussian-like" basis set that is confined to the trun- 
cated space. We use the basis sets {gi(x)} and {0j(x)} 
to construct a new basis set, {cji(x)}: 



N 



9m{x) = ^ d n (x)g m (x n ) 
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for m = 1, N. The new basis set is in some sense, the 
Gaussian functions with periodic boundary conditions. 
We can write eq.Q in matrix notation as: G = 0G 
where = gj{xi) By taking the width parameter a — 
we can guarantee that the pvN functions have no 
linear dependence and that the matrix G is invertible, 
that is GG^ 1 — 9. The invertibility of G implies that 
both bases span the same space. 

The representation of in the pvN basis set is given 
by: 
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To find the coefficients a m we first define the overlap 
matrix, S: 



{9i\9j) = I 9*{x)9j(x)dx 
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FIG. 1: N = 9 coordinate grid points and N = 9 vN unit cells 
span the same area in phase space, S = Nh. The vN basis functions 
are Gaussians located at the center of each unit cell. 



contributions from remote basis functions, but we can 
take advantage of the vanishing overlaps by defining a 
bi-orthogonal von Neumann basis (bvN) {bi(x)}. 
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or in matrix notation: B = GS 1 . Inserting eq[9] into 
eqJSl can be written as 
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S = G 1 G. 
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Using the completeness relationship for non-orthogonal 
bases, can be expressed as 
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Comparing with eq.© we find that a m = 
Y,n=i( S ^ 1 )mn(9n\i}) and (g t \^) =T, n =l9i( x ^( x n)- 

Although the periodic von Neumann (pvN) and the 
Fourier methods span the same rectangle in phase space, 
the localized nature of the basis functions in the pvN 
method can lead to significant advantages. In particu- 
lar, if \tp) has an irregular phase space shape we may 
expect that some of the pvN basis functions will ful- 
fill the relation: (gj\ijj) = 0, j = 1,...,M. Due to the 
non-orthogality of the basis we cannot simply eliminate 
the states cjj, since the coefficients of cjj may include 



By assumption, M of the coefficients are zero, hence in 
order to represent \tp) in the bvN basis set we need only 
N' = N~M basis functions. Note that the bvN and pvN 
are bi-orthogonal bases, meaning that each set taken by 
itself is non-orthogonal but they are orthogonal to each 
other. This can be shown easily by: 

JV 

(9i\bj) = ^9i{Xn)fj{x n ) 

71=1 
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For many practical applications the full knowledge of the 
basis wavefunctions is unnecessary: we need only the 
value of the basis functions at the sampling points. For 
example the evaluation of Hamiltonian matrix elements 
can be performed explicitly by: 
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and similarly: 
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where H FGH = V FGH + T FGH and the potential and the 
kinetic matrix are given by: V? GH « V(xi)8ij and 



T; 



FGH 



+ if i=j 
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[16J. The eigenvalue problem in a non-orthogonal basis 
set becomes HU = slIE; in the pvN basis set s is given 
by eq. and in the bvN basis set s is given by: 



B*B = S-i&GS- 1 = S- 1 . 



(15) 



Diagonalization should give accurate results for all wave- 
functions localized to the classically allowed region of the 
rectangle. 

As a numerical test of the pvN basis set we studied 
the standard example of the harmonic oscillator V(x) = 

2 2 

m " 2 x in units such that m = n, = oj = 1. We calcu- 
lated the first 8 eigenenergies using 16 pvN compared 
with 16 conventional Gaussian basis functions. In the 
Gaussian basis set the Hamiltonian and the overlap ma- 
trices were calculated analytically as: ffy = (gi\H\gj) = 
JZo9*( x ){-7& + V(x)]g j (x)dx and 5 y = ( gi \gj) = 
J '_ g*(x)gj(x)dx. The 16 sampling points were taken 
from -5 to 5-S x and the width parameter was a = 0.5. 
The results, shown in Fig. [3J show the superiority of the 
pvN basis set compared to the standard Gaussian basis 
set. In fact, the results obtained with the pvN basis set 
are exactly as accurate as in the Fourier grid method. 
The kinetic energy spectra in Fig. [5]reveal the deficiency 
of the conventional Gaussian scheme. 

In the bvN basis set we are able to remove some of the 
basis functions and construct lower dimensional H hvN 
and 5' bvN matrices without losing accuracy. In order to 
test this claim, we calculated numerically the eigenener- 
gies of the Morse oscillator V(x) = D(l — e~@ x ) 2 and the 
Coulomb potential V(x) = by using both the FGH 
and fVN basis sets. The Morse parameters were taken 
to be D = 12, m = 6, f3 = 0.5 and H = 1. For FGH, 
100 grid points between [—1.6, 20.1] were required in or- 



o °o Q o 

CD 

—5 A 
CD ▲ 
O A 



-10^ 



10 . 




3 6 

n 



FIG. 2: (a) Error in the lowest 8 eigenvalues of the harmonic 
oscillator. vN(O), pvN(A). (b) Kinetic energy spectra using 16 
basis functions. vN(dashed), pvN(solid). 



der to get 4 digits of accuracy in energy for all 24 bound 
states. By using the bvN basis functions (constructed 
from 10x10 vN functions with a = 0.5) we obtain the 
same 4 digit accuracy with only 48 basis functions. This 
is demonstrated graphically in Fig. [3] (a). The figure 
shows the phase space representation of 100 evenly grid 
points. Although it requires 100 pvN basis functions to 
span this area in phase space, due to the flexibility of the 
bvN basis set we can suffice with just the basis functions 
in the classically allowed region (white squares). 
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FIG. 3: Phase space area spanned in the bvN method (white) and 
in the pvN (or FGH) method (full rectangle) for Morse (a) and 
Coulomb (b). 



In the Coulomb potential the efficiency of the bvN ba- 
sis set is even higher. The rectangular shape in phase 
space of the FGH is a very wasteful representation for 
the cross-like shape of the eigenstates. The Coulomb pa- 
rameters were taken to be Q = 1, m = 1 and h = 0.5. 
For FGH 1599 grid points between [—50.2, 50.1] were re- 
quired in order to get 4 digits of accuracy in energy for the 
first 9 excited states. By using the bvN basis functions 
(constructed from 39x41 vN functions with a = 0.2367) 
only in the classically allowed energy shell we obtain the 
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same accuracy with only 189 basis functions. Fig. [3] 
(b) illustrates the efficiency of the bvN basis by showing 
the phase space area in FGH and bvN bases. Figure @] 
shows the error in the eigenenergies by using 190 func- 
tions in the FGH and in the bvN method respectively. 
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FIG. 5: Basis efficiency for the Morse potential as function of ft. 
For pvN (or FGH) (dashed) and bvN (solid) 
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FIG. 4: The deviation between the calculated and the exact 
eigenenergies for the lowest 20 eigenstates of the Coulomb potential 
using 190 basis functions for the bvN (a) and pvN (or FGH)(0) 



The ability to localize a bvN function at a specific point 
in phase space results in the remarkable concept of 1 ba- 
sis function per 1 eigenstate. This means that in order 
to calculate N eigenenergies we need only N basis func- 
tions. Obviously, such one per one efficiency, if reachable, 
will be the ideal efficiency for any basis set. In order to 
test the ability of the bvN method to reach the ideal 
efficiency we examined the Morse potential and looked 
for the smallest number of bases that will provide ex- 
act values of the energies (12 digits of accuracy) for all 
the eigenstates in the energy shell E = 11.25. The bvN 
method indeed, tends to the ideal efficiency in the clas- 
sical limit where h — > (Fig. [5]). This remarkable result 
is unique for the bvN method. For example, the effi- 
ciency of the Fourier method has an upper bound, which 
can be calculated analytically as 77 = 1.6959 (the ratio 
of the classically allowed phase space to the area of the 
rectangle). 

In summary, we introduced the "Gaussian-like" pvN 
basis set and showed the equivalence between the pvN 
and the Fourier method. Several TISE problems were 
tested. We showed that the pvN method is as accurate 
as the Fourier method; however, due to the localization of 
the basis functions, we can construct a flexible basis set, 
bvN, which is much more efficient then the FGH. We also 
have shown that in the classical limit the number of basis 
functions is equal to the number of the exact eigenvalues. 
This remarkable "1 per 1" quality, is of course, an optimal 
result for any basis set, and as far as we are aware has 
never been achieved before. 

As discussed above, the use of primitive Gaussian 
functions on the vN lattice yields disappointing results 



12j, [l3j. The success of our method is due to two key 
points: 1. By defining the basis set as a linear combina- 
tion of the {8i(x)} we ensure that the basis indeed spans 
a band limited region with finite support. 2. By using 
the bvN basis set the coefficients become locallized. In 
the conventional Gaussian basis set the delocalized coef- 
ficients arise from the non-locality of S^ 1 . 

In this paper we focused on model 1-d TISE problems. 
However, Gaussians are used as basis functions in a va- 
riety of methods for solving the TISE and TDSE[3-[2l|. 
All these methods have suffered from numerical difficul- 
ties resulting from the non-locality of S* -1 and we be- 
lieve that significant improvement is possible using the 
pvN and the bvN basis sets. Finally, the pvN method is 
not limited to quantum mechanical problems. There is 
a large literature on the use of a Gaussian basis set in 
signal processing where it goes by the name of the Ga- 
bor transform [22h26| . The Gabor transform is known to 
have problems with stability, which can be traced to the 
"no-go" theorem of Balian and Low[27l 28 1 — the state- 
ment that localized basis sets are incompatible with or- 
thogonal basis sets. While the pvN and bvN bases do 
not violate the no-go theorem, they seem to effectively 
circumvent it. We therefore believe they can have a sig- 
nificant impact on signal processing in general. 

This work was supported by the Israel Science Founda- 
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